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Qh! Abstract 

(D 

i-^ ' We determine the bottom MS quark mass rrlh and the quark mass in the po- 

tential subtraction scheme from moments of the hh production cross section and 
KA \ from the mass of the Upsilon IS state at next-to-next-to-leading order in a reor- 

^H ' ganized perturbative expansion that sums Coulomb exchange to all orders. We 

"^ ' find mb{mb) = (4.25 ± 0.08) GeV and mfe,ps(2GeV) = (4.59 ± 0.08) GeV for the 

potential-subtracted mass at the scale 2 GeV, adopting a conservative error esti- 
mate. 



Introduction. Accurate determinations of the bottom quark mass in perturbative 
QCD usually rely on properties of the spectrum of Upsilon mesons and bb production near 
threshold. Since already for the T(1S) state the momentum scale p = mbCpasip)/'^ ~ 
1.25 GeV and energy scale E = mb{CFas{p)Y/4: < 0.3 GeV are too small for perturba- 
tion theory to be expected to to work, one considers averages over the bb production 
cross section through a virtual photon, including the T resonances^ [Q: 

Mj(lOGeV)- ^ -^^^.^i<l\.o = /^^-(^)- (1) 

In this case the characteristic momentum and energy scales are replaced by p ~ Ivn^yj ^fn 
and E ~ mb/n, respectively. The requirement of perturbativity puts an upper limit on 
the admissible values of n. On the other hand, only the resonance contribution to the 
sum rule is experimentally well known, and n needs to be taken large enough to reduce 
the error from the bb continuum. When n ^ 1, the perturbative expansion of M„ in the 
strong coupling as breaks down, because there exist terms of the form {as\/n)^ in any 
order of perturbation theory. This suggests a summation of the perturbative expansion 
to all orders in which a^-y/n is treated as order 1 [^. 

In this letter we analyse the sum rule (|l|) at next-to-next-to-leading order (NNLO) in 
this resummed perturbative expansion. [A preliminary analysis was presented in Ref. |^.] 
The resummed perturbative hb cross section is computed at NNLO using recent 2-loop 
results on the Coulomb potential |Q and the 7*66 vertex [^ H, and non-relativistic 
effective field theory in dimensional regularization as described in |^, ^, ^. Rather than 
determining the b quark pole mass, as has usually been done, we apply the potential 
subtraction (PS) scheme and determine the PS mass |]10[ from the sum rule. We expect 



perturbative corrections in this and related schemes to be smaller than in the on-shell 



scheme [[L0| |Tl|. We then convert the extracted PS mass to the MS mass, thus by- 



passing the infrared sensitivity problem of the on-shell scheme |]T2[, and yet implementing 
the resummation necessary in the non-relativistic kinematics enforced by taking large 
moments. Other NNLO analyses of the sum rule have already appeared [|1^, [1^, |15|, ^ . 



Nevertheless, we think that an independent analysis, together with a critical discussion 
of the quark mass error, is still useful. We also perform a complementary analysis and 
determine the b quark mass directly from the mass of the T(1S) state. This has been 
done previously in a NNLO analysis presented in []T7|, which, however, concentrated on 
the b quark pole mass, as did |]13|, |14|. 



Experimental moments. We first evaluate the integrals (0) by expressing the cross 
section in terms of the six T resonances and the open bb continuum. The masses and 



leptonic widths of the resonances are taken from [18|. Very little information exists 



on the bb continuum above s = (10.56 GeV)^ |l^. We parametrize the continuum by 



setting Rli = 0.4 ±0.2. With this crude parametrization the experimental error on the 
determination of m^ is 110 (50, 30, 15) MeVforn = 4 (6, 8, 12), and small compared to the 
theoretical error for interesting moments with n ~ 8-12. Some experimental moments 



^ The cross section is normalized such that Rhiis) = Nccl = 1/3 in the ultra-relativistic hmit, where 
Cfc is the bottom quark electric charge. Ri,i{s) is related to the vector current two-point function n(g^) 
in the usual way. The normalization factor on the left-hand side of (yj) is inserted for convenience. 



n 


4 


8 


10 


12 


16 


M„ 


0.231 ± 0.044 


0.235 ± 0.023 


0.264 ± 0.021 


0.308 ± 0.021 


0.443 ± 0.025 



Table 1: The experimental moments. 



are shown in Table |I]. For n 
from the T(1S) resonance. 



8-12 about 70%-85% of the experimental moment comes 



Theoretical moments. The theoretical moments are computed by first matching 
QCD to non-relativistic QCD. In a second step this theory is matched to a non-local 
Schrodinger field theory, in which hh pairs propagate through the Coulomb Green func- 
tion. We then solve the Schrodinger equation to NNLO. We refer to 0, ||, |^ for some 
details of the method; further useful information can be found in |l^, |14|, |15|. The result 
for the hh cross section to NNLO, still in the on-shell scheme, is expressed as 



R, 



■bb\ 



as{fo{XJ) + asfi{X,l)+alUX,l)} 



(2) 



where / = \n{-AmbE / /j,^) , X = CFa,/(2 {-E/mbY/^), Cp = 4/3, E = ^/^-2mb and rub 
is the h pole mass. The functions fi{X,l) contain bound-state poles that correspond to 
the T(nS) resonances. We obtained these functions analytically. After integrating nu- 
merically over s according to (|1|), these functions sum all terms of the form a°'^'^ (c^sV^)'^ 
to all orders. Writing Rbb{s) in the form of (H) implies that we expand the bound-state 
pole (5-functions around the leading-order pole position. Expanding the bound-state pole 
^-functions rather than leaving them unexpanded is motivated by the fact that the sum 
rule relies on global duality. Using dispersion relations, the moments can be expressed 
in terms of derivatives of the vacuum polarization as indicated in (|l]), which makes no 
reference to individual resonances. Computing these derivatives in resummed perturba- 
tion theory to NNLO implies that we expand the resonance (5-functions in the expression 
fori?,5(s).g 

Before integrating over s, we convert the expression for Rbb{s) from the on-shell to 
the potential subtraction scheme. The pole mass is eliminated using the relation |]10 



rrib = mb^psinf) - - 



"^b,ps(/i/) + 



|g|<M/ 
Cpa 



(2^ 



Viq) 



TT 



7^/ 



as 
Air 



6mi + 



a 



^ ] 6m2 + . . 
Air/ 



(3) 



where V{q) is the Coulomb potential in momentum space. Explicit expressions for 6mi^2 
can be found in [^. Note that rrib — rrib^psif^f) is proportional to a subtraction scale ///, 
which should not exceed the characteristic scale 2mb/y/n of the moments M„. We insert 
(§) into @ and expand the small correction terms involving 6mi^2- However, the term 

^ For very large n one should keep the ^-functions unexpanded, because the effective smearing interval 
in s becomes smaller than the perturbative correction to the bound-state pole position. But for such 
large n one has to rely on local duality and the sum rules suffers from non-perturbative uncertainties as 
we discuss further below. 



A = Cpasfif/iT is not expanded when rrib is replaced in E, X or /, because A counts as 
being of the same order as i^^ = y/s — 2mb. The result is an expression of the same form 
as (H), but with mb^ps(/i/) as input parameter. As mentioned in the introduction we 
expect the expansion (0) in this new variable to be more convergent, and hence the PS 
mass can be determined with smaller error than the pole mass. The PS mass depends 
on Hf and we choose fif = 2 GeV as our default. The PS masses for different jU/ are 
connected by a renormalization group equation that follows directly from the definition 

(I)- 

The dominant theoretical uncertainty arises from the residual dependence of the 
theoretical moments on the renormalization scale /z of the strong coupling a.^ = dsif^) in 
the MS scheme. We now discuss the choice and variation of this scale and the choice of 
moments n that go into our analysis. 

As indicated in (^ explicit logarithms of /i always come as ln{—4:mhE / fi'^) . When 
n ^ 1, the integral (|I]) falls exponentially as exp{—nE/mh), so that the characteristic 
energy scale is i? ~ mf,/n. This determines the parametric form of the renormalization 
scale to be /i„ = 2mb/y/n. This, however, is not strictly true, because the moments also 
contain parts in which gluons carry momentum of order rrih and momentum of order 
iTih/n. In the renormalization-group- improved treatment (see 0, P) the hard scale mj, 
enters as the starting point iih ~ '>Tib of the renormalization group evolution of the Wilson 
coefficient functions. The dependence on /x/j is negligible, of the order of ±10 MeV on 
the output for m^, compared to the dependence on the scale /i„, which determines the 
endpoint of the renormalization-group evolution. It is therefore not considered further. 
Gluons with three-momentum of order r/ife/n enter only at order a^. This leaves us with 
the scale /i„ above and we adopt /i„ = 2mh/ \/n as the most 'natural' scale. 

One may object that the form of the logarithm gives the natural scale only paramet- 
rically, but that the scale is arbitrary up to a multiplicative factor, since the physical 
scale in the MS scheme corresponds to a different scale in another scheme, for instance 
MS. We can address this question by searching for constants that appear systematically 
in conjunction with the logarithm ln(— 4?7ib£'//i^). While this is complicated for the 
full cross section, it is easily done for the bound state energies that correspond to the 
T(nS) resonances. We find that for the /cth energy level the analogous logarithm always 
appears in the combination \n{mi,C pds / {kfj)) — Si{k), where Si{k) = J2m=i V""^- Since 
Si{k) > 0, this suggests - if anything -, that the physical scale is even smaller than what 
we inferred from the logarithm alone. 

The useful moments are restricted from below by the uncertainty in the experimental 
value of the moment. If we aim at an error of about 50 MeV in the determination of 
the b quark mass, we need n > 6. There is also a technical restriction, which could be 
overcome. Our expression for M„ sums all terms of the form a^'^'^(a;sV^)'^, but it does 
not make use of the exact fixed-order coefficients at order a]'"^, which are known [ ^ , 
because terms of relative order n~'^/^ or smaller are dropped. This could be compensated 
for by matching the resummed result and the fixed-order result. However, we find that 
for n > Q this matching correction is small, as can be seen from Table 0. 

It is advantageous to take large moments, because large moments are more sensitive 
to rrib, while the experimental error does not increase, see Table |I|. An upper limit 
arises, because the characteristic scales must remain perturbative. As concerns /i„, the 





LO 


NLO 


NNLO 


n 


Res. 


FOl 


F02 


Res. 


FOl 


F02 


Res. 


FOl 


F02 


6 


0.134 [0.228] 


0.030 


0.028 


0.127 [0.143] 


0.062 


0.065 


0.201 [0.222] 


0.094 


0.097 


10 


0.122 [0.268] 


0.016 


0.016 


0.143 [0.166] 


0.041 


0.042 


0.266 [0.322] 


0.073 


0.075 


14 


0.139 [0.379] 


0.011 


0.011 


0.190 [0.218] 


0.032 


0.033 


0.400 [0.529] 


0.065 


0.067 



Table 2: Selected moments (recall the normalization of the moments in (|lj)) in the on-shell 
and PS scheme (square brackets) for as{Mz) = 0.118 and rrn, = 4.95 GeV (on-shell scheme) 
and mb^ps(2GeV) = 4.55 GeV (PS scheme). The renormalization scale is taken to be ^„ = 
2M/\/n, where M = m\, or M = TTifo^pg (2 GeV) . 'Res.' refers to the resummed cross section 
with LO, NLO, NNLO in the sense of @). 'FOl' refers to the fixed-order result without 
resummation, including terms of order a^'^'^ (LO, NLO, NNLO). 'F02' refers to the fixed-order 
result, dropping terms with relative suppression n~^' ^ or more. The effect of resummation is 
given by the difference of 'Res.' and 'F02'. The correction due to matching to the fixed-order 
result is roughly given by 'FOl' minus 'F02'. 



requirement /i„ > Aqcd does not seem to pose a serious restriction. In practice, we find 
that the theoretical prediction becomes unstable already when /i is smaller than 1.5- 
2.0 GeV; requiring /i to be larger than this is restrictive, if we also allow for a variation 
of /i about /i„. A more serious constraint arises from the scale mb/n, which enters the 
NNLO calculation implicitly. At NNNLO there is a contribution to the moment that 
scales as as{iin)'^(^s{.''^h/n)- When m^/n ~ Aqcd, we should count as{mh/n) as order 1. 
In this case, we have an uncontrolled non-perturbative contribution to the moments that 
is formally of NNLO.0 We therefore require n < 10. In the literature larger moments are 
often used. The justification for this is that the gluon condensate contribution to the 
moments, which represents the leading non-perturbative power correction, is small even 
for moments much larger than 10. However, the operator product expansion in local 
operators is itself only valid when mf,/n > Aqcd and so the estimate is not rigorous. It 
may, however, indicate that ultrasoft contributions from the scale rrih/n are smaller than 
what we would estimate on parametric grounds. 

For reference we give some selected moments in Table in the on-shell and PS scheme. 
The table also quantifies the importance of resumming {as\/n)^ corrections and the error 
incurred by not including the exact fixed-order coefficients at order a^'^. Resummation 
is crucial even for n = 6, because the contribution from the bound-state poles, which 
does not exist in the NNLO fixed-order approximation, is large. On the other hand, 
already for n = 6 the fixed-order moments (FOl) are well approximated by the leading 
three terms in their large-y^ expansion (F02). 

Numerical analysis. For a given n the theoretical moments are functions of r7i5Ps(yUj), 
which we would like to determine; the strong coupling as, for which we use as{mz) = 

■^It is worth noting that the scales /i„ and mb/n do not really approach as n — > oo, but freeze at 
values of order mtas and rriba'^, respectively. While this is of interest for a very heavy quark, it is of 
little practical relevance to b quarks. 
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Figure 1: The value of mb^ps(2GeV) obtained froui the 10th moment as a function of the 
renormahzation scale in NLO and NNLO and for as{mz) = 0.118. The dark region specifies 
the variation due to the experimental error on the moment. The middle line marks the scale fin, 
the two outer lines determine the scale variation from which the theoretical error is computed. 




0.118 ± 0.003 together with 3-loop evolution; and the renormalization scale /i, which is 
our (rough) handle to estimate the uncertainty due to the NNLO approximation of the 
moment calculation. We first compute, for given n and /i (and agimz) = 0.118), the val- 
ues of mf, ps(2 GeV), for which the theoretical moment lies within the experimental range. 
For n = 10 the result is shown in Fig. ^ It is evident that the resulting m6_ps(2GeV) 
varies significantly as function of the scale /x at which the sum rule is evaluated. Further- 
more, there is no overlap between the range of masses that is obtained from the NNLO 
and the NLO sum rule for any reasonable range of /i. The same conclusion is obtained 
for n = 6 or 8. We also determined ?Tife^ps(2 GeV) from n = 6,8, 10 simultaneously by 
minimizing a x^ with equal weights. The value we obtain from this procedure differs 
by no more than 30 MeV from that obtained from single moments, when /i is varied 
between 1.5 and 9.5 GeV, reflecting that the /^-dependence of the theoretical moments is 
completely correlated. The same analysis in the on-shell scheme results in an identical 
qualitative picture; however, the scale-dependence is even larger in the on-shell scheme. 
As explained above, we take /i„ = 2mfe,ps(Ai/)/\/n as our default choice of scale. We 
would then follow common practice and estimate a theoretical error by varying the scale 
between one half and twice this value. But from Fig. [^ we observe that the theoretical 
prediction becomes unstable (compare the behaviour of the NLO and NNLO results) 
for scales below 2 GeV and one may argue that varying the scale into this region does 
not provide a reliable error estimate. We therefore compute the theoretical error from 
a variation between 2 GeV and 2/i„. It is clear that the error so estimated is rather 
sensitive to the lower scale cut-off. Taking n = 10 and adding the error from as and the 
experimental moments, we obtain 



mfe,ps(2 GeV) = (4.60 ± 0.10 (scale) ± 0.03 (a^) ± 0.02 (exp.)) GeV. 



(4) 



If the scale is varied down to 2.5 (1-5) GeV, the scale error decreases (increases) to 
± 65 (160) MeV. In comparison, a NLO analysis of the sum rule would return the central 



value mf,^ps(2GeV) = 4.44 GeV with a smaller scale uncertainty (see Fig. H). The large 
difference with the NNLO result casts doubt on the convergence of successive pertur- 
bative approximations. The origin of this difference and the origin of the large scale 
dependence will become clear below. 

The PS mass is a useful parameter (replacing the pole mass) for short- distance ob- 
servables involving b quarks close to their mass shell. For high energy processes, we 
would like to convert the PS mass to the MS definition. Call nib the MS b quark mass at 
the renormalization scale fni, and kr the coefficient at order {as{rnb)/{4:'^)Y that relates 
the pole mass to m^. From (13) we obtain the relation 



/ r ^^f 

mb.PsifJ'f) = rnb\l + Y, ^r - 4:CFSmr-i{fif) ^ 

\ r=l L ^b 






where we defined 6mo = 1. An N'^LO analysis of the sum rule determines the PS mass 
with a parametric accuracy of order rub^ps a^^^- This is most easily seen by noting that 
an N'^LO calculation of the resummed bb cross section determines the (perturbative) 
T(nS) masses to order mbpsa^"^^. To determine the MS mass with the same parametric 
accuracy implies that one should use (|^) at order aj"*"^. At present (^) is known only to 
third order, combining the result of [^] for k^ and the one for 6m2{fif) from pO| . 



To obtain an order-of-magnitude estimate of the missing fourth-order term, we es- 
timate the coefficient in the 'large- /3o' approximation P2|, EH]. This gives k^ = 339071 



and (5m3(2GeV) = 125026 for Uf = 4 and Mb = 4.26GeV.g Although the individual 
coefficients are large, there are large cancellations in the combination that enters (|^), 
which reflect the infrared cancellation that motivated the introduction of the potential 
subtraction [0. With these numbers, given m;,ps(2 GeV), we estimate that the a^ term 
reduces m^ by 10 MeV. (For comparison, the a^ term provides a 35 MeV reduction.) We 
therefore assume an additional (10 ± 20) MeV correction in the relation between mb and 
m6,Ps(2GeV) beyond the third-order formula. This results in the MS massQ 

nib = (4.26 ± 0.09 (scale) ± 0.01 (a,) ± 0.02 (conv.) ± 0.02 (exp.)) GeV. (6) 

The dependence on as nearly cancels out and 'conv.' refers to the conversion from the 
PS to the MS scheme just discussed. We have repeated the analysis with /ij = 1 GeV 
and 3 GeV for the subtraction scale of the PS mass. Converting to nib, we find agreement 
with d) within 20 MeV. 

Origin of the large scale dependence. The scale uncertainty in (|) is only about 30% 
smaller than the uncertainty we would have found in the on-shell scheme. To understand 
the origin of this marginal improvement, we consider a truncated sum rule, in which both 
experimental and theoretical moments are given only in terms of the first T resonance. 

''For comparison, note that the 'large-/?o' approximation for 6m2{2GeV) results in 2140.36 rather 
than the exact value of 1870.54. For k^ one obtains 6526.91 rather than the 'exact' value 6144(128). 
(The brackets specify the error on the 'exact' result, see |2^.) nj refers to the number of light-quark 
flavours. 

^ We compute rnb by solving (g|) exactly for a given PS mass, rather then inverting (J^) perturbatively 
to order a*. The second procedure would result in a central value that is negligibly different by 4 MeV 
from the one given. 



This is actually not a bad approximation to the full sum rule and allows us to discuss 
the origin of scale dependence in a transparent form. In this approximation, performing 
in addition a non-relativistic approximation to the s-integration measure in (0), we can 
write the sum rule as 



r*i._ r . . M^^r., - Mi!" 



l+l- _ 

exp ^^P 



v 



{2n + 1) 



^T( lS) ^"Y(1S) 



(7) 



M^lis) = 2mb - - TUbai 1 + [3.590 - 1.326 /] «, + [19.52 - 6.033 / + 1.319 r 



a„ ., 



where rJ+;_ and M^'/^g) are the leptonic width and mass of the T(1S) state computed to 
NNLO. In the on-shell scheme, the series expansions for the two quantities readQ 

4 / 

rth _ o™ 2 / 1 , ro r^nn i ooc ;1 „ , n n co a noo i , i 01 n r2i „,2 

9 

= 2mb--mbal{l + im + l.7Q + ...\ (8) 

32 / 

rfi;_ = — elal^m^al 1 + [3.605 - 1.989 / - 5.607 k] a, 

+ [24.14 - 7.898 / + 2.639 f - 0.3793 In k - 20.21^; + 11.16 /k + 7.260 k^] a 

32 / \ 

= ^ eWe^nibal [l - 0.11 + 1.23 + . . . j (9) 

where a^ = ^^(/u), / = ln(16m^a^/(9/i^)) and n = as{m)/as{fi) and the second line is 
given for /i such that / = (for rrih = 5 GeV), in which case as(/i) = 0.30. Neither of the 
series is converging and, because of the large power of a^ of the NNLO term, the scale 
dependence is huge at small scales, where a^ varies fast. This is seen from Fig. |^, which 
shows the scale dependence of the left-hand side and right-hand side of (|^ separately. 

Going from the on-shell to the PS scheme improves the convergence and scale de- 
pendence of the predicted T(1S) mass as seen from the solid line in Fig. |^, but has little 
effect on the series expansion (^ for the width. Hence the large scale uncertainty in (|) 
can be traced to the poor control over the perturbative expansion for the leptonic width 
that controls the over-all normalization of the theoretical moments. 

Constraints from Mx(is)- The fact that the mass determination from the sum rule is 
dependent on the theoretical prediction of the leptonic width, and limited in accuracy 
for this reason, suggests that we consider determining rrih directly from the T spectrum. 
Non-perturbative corrections to the T masses grow rapidly for higher radial excitations 
and preclude using any state other than the ground state. The problem with this method 
is that even for the T(1S) state it is difficult to estimate the non-perturbative correction 
reliably. The perturbative expression for the T(1S) mass in the on-shell scheme is given 
by (ID above. 

^The mass and leptonic width have been obtained to NNLO in ||l^, ll^ and fl^, respectively. Our 
analytic expressions for an arbitrary T(nS) state coincide with those previous results, provided we 
neglect the renormalization-group improvement for the leptonic width. At present only the logarithms 
from the renormalization of the external current are taken into account. 
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Figure 2: The left-hand side (Ihs, short-dashed line) and right-hand side (rhs, dash-dotted 
and solid lines) of (0) as a function of the scale /i for n = 10, m^ = 5 GeV (on-shell scheme, 
dash-dotted) and rn,ft.ps(2GeV) = 4.6 GeV (PS scheme, solid). The figure clearly shows the 
reduction of the scale dependence in the PS scheme for the predicted T(1S) mass. The scale 
dependence of the width (short-dashed line) is identical in the on-shell and PS schemes. 



If mb{CF(ysY/4: 3> Aqcd, the non-perturbative correction to the T(1S) mass can be 
computed in terms of vacuum condensates of local operators. The leading contribution 

is Pl 

_ 624 _ {asGG) 



5M. 



T(1S) 



425™'(mbCi.a,)4^ 



(10) 



where {agGG) ~ 0.05 GeV^ is the gluon condensate. The actual magnitude of 5M^,^^\ 
is rather uncertain. If we choose the 'natural' scale /z = mhGpasifi) ~ 2 GeV, we obtain 
5M^?-^g^ ^ 90 MeV. However, as noted earlier, the logarithm that determines this scale 
appears together with constants that tend to make the effective scale lower. Furthermore, 
it may be argued that the coupling as should be taken as the coefficient of Cp/r in the 
Coulomb potential rather than in the MS scheme. This coupling is larger than the MS 
coupling. Both effects can decrease (5M^^gx substantially. 

Since the inequality mb{GFasY/4: 3> Aqcd that justifies the operator product expan- 
sion (OPE) does not hold, we should consider the subsequent term in the OPE to judge 
whether the expansion converges. Using the result of [^, we find that the contribution 



from dimension-6 operators could be anything between a fraction of and twice — (5M^^gN, 
where the large uncertainty stems from the poorly known dimension-6 condensates and 
the ambiguity in the value of ag!^ This puts the convergence of the OPE in question. 
We therefore consider (|TD|) as an order-of-magnitude estimate of the non-perturbative 
correction and treat it as a theoretical error rather than adding it to (|]). In our opinion, 
assigning an error of ±70 MeV to nib from this source is conservative. 

We proceed to determine mf,,ps{fif) from Mr{is) using (||) converted to the PS scheme. 
This renders the series (§) convergent and leads to a very small scale uncertainty of the 
extracted value of mb^ps(2GeV), as anticipated from the solid curve in Fig. ^ Varying 

''^Ref. 12^ concludes that the OPE appears to be convergent, because the minimal value of the strong 
coupling in the denominator of (|l^ assumed there is larger than the minimal value allowed in our 
estimate. 



/i from 1.25 to 4GeV, we obtain 

mb,ps(2 GeV) = (4.58 ± 0.04 (scale) ± 0.01 (a,) ± 0.07 (non-pert.)) GeV, (11) 

which is consistent with (^). In this case a NLO analysis would return the central value 
m6_ps(2 GeV) = 4.53 GeV, which suggests that the corresponding small value in the sum 
rule analysis is an anomaly related to the behaviour of the series for the leptonic width. 
From (1TT|) we obtain the MS mass 

Wb = (4.24 ± 0.04 (scale, /i/) ± 0.02 (a,) ± 0.02 (conv.) ± 0.07 (non-pert.)) GeV. (12) 

The central value varies by only about 10 MeV, when fif is varied between 1 and 3 GeV. 
Contrary to the sum rule determination, the error is dominated by the non-perturbative 
contribution to the T(1S) mass. This leaves room for improving upon the error, if some 
quantitative insight into the non-perturbative contribution could be obtained. 

Comparison with previous results. We compare the bottom quark mass obtained in 
this work with the results of earlier NNLO analyses of the sum rule and the T(1S) mass. 
Our comments will be restricted to those analyses that quote a result for the MS quark 
mass 



13| , p!4| , p!5| , p!6| , p!7| , |26| . With the exception of [|13], which obtains m;, from Mx(is) 
and, therefore, should be compared with (p!2D, all other NNLO analyses use the sum rule 
(HI) and should be compared with (|^). 

The value given in (|T2p is significantly smaller than m^ = (4.44 ±0.04) GeV, obtained 



by 113 • This difference is explained by the fact that [^ first uses the on-shell scheme to 
extract the pole mass and then uses the 2-loop truncation of (|^) [with fif = 0] to obtain 
nib- However, contrary to the PS scheme with fij = 2 GeV, the 3-loop and 4-loop terms 
are large in the on-shell scheme; at least the 3-loop termQ has to be included when the pole 
mass is determined from the NNLO formula for the T(1S) mass. Estimating the terms 



missing in [JT3 in the large-/3o limit, and subtracting them from nib = (4.44 ±0.04) GeV, 



we find that the result of [T^ becomes (roughly) consistent with ours. The error estimate 



of [|T^ is, however, less conservative than ours. 

A related difficulty concerns the comparison with the value m^ = (4.21 ± 0.11) GeV 
quoted in [O]. While apparently consistent with the one obtained in this work, it is 



obtained via a 2-loop relation from the b quark pole mass, which in turn is determined 



from the sum rule. If we add the 3-loop and/or 4-loop term, the result of [13| would be 



about 200 MeV lower than ours. This difference is a reflection of the fact that the pole 



mass quoted in [|13| is roughly 200 MeV lower than the one we would have obtained had 
we chosen to determine it. This difference in turn can be traced to the use of a high 
renormalization scale for the evaluation of the sum rule, cf. Fig. 0. In our opinion, the 
choice of such a high scale is not well motivated. We also think that it is mandatory to use 
intermediate mass definitions such as the PS mass to determine nib reliably. Otherwise 
large perturbative coefficients make it impossible to disentangle true theoretical errors 
from correlated and spurious ones caused by those large perturbative coefficients. 



^Since the large infrared contribution in the 4-loop term cancels against an NNNLO contribution to 
the T(1S) mass, it can be argued that only the 3-loop term is to be used. This is different from the 
PS scheme, where no systematically large coefficients appear. Compare (^) [with /x/ ~ nibas] and the 
discussion in [|lO| , |27| regarding combining different powers of as to make infrared cancellations manifest. 



The most recent analysis by the authors of [^] determines the MS mass directly 



from the sum rule and gives mj, = 4.19 ± 0.06 GeV from moments with n = 7-15. No 
resummation is performed, because it is assumed, incorrectly, that this is unnecessary in 
the MS scheme. However, for high moments the Coulomb interaction must be treated 
non-perturbatively, and a resummation has to be done, irrespective of the mass renor- 
malization convention. The MS scheme actually makes the expansion worse, because the 
expansion contains terms of order (a^n)'^ in addition to (asy/n)^. To avoid such terms, 
one has to use an intermediate convention, such as the potential subtraction scheme, 
and then relate this convention to the MS scheme in a second step. Because of this 
theoretical shortcoming, the result of [^ cannot be compared with (^). 



The other papers quoted above perform a NNLO resummation as in this work. Dif- 
ferences arise either in the representation of the NNLO-resummed moments or in the 
analysis and error evaluation strategy. Both [|I^ (MY) and |T^ (Hoang) also use inter- 



mediate mass subtractions, different from the PS scheme, but conceptually similar to it, 
before converting these intermediate masses to the MS scheme.0 

The differences in the theoretical representation of the moments are the following: 

(a) MY and Hoang use a factorization scheme different from dimensional regulariza- 
tion. Since the final result is physical, this is a technical difference that should 
bear no consequences on the final result. 

(b) Hoang applies the non-relativistic approximation also to the s-integration in (|l|), 
while MY and the present work obtain the resummed cross section analytically 
and then integrate it numerically according to (|1|) or after an equivalent contour 
deformation into the complex s-plane. The difference is negligible. 

(c) MY and Hoang have taken the short- distance coefficient as an over-all factor, while 
we have multiplied it out to NNLO. Keeping it as an over-all factor is problematic, 
because this results in a spurious factorization scheme and scale dependence, which 
is not small as can be seen from Table 3 of |T6|. Since both short-distance and 



long-distance contributions are computed perturbatively, the factorization scale is a 
purely technical construct and no dependence on it should be left in the result. One 
motivation for writing the short- distance coefficient as an over-all factor is that the 
scales in the coupling constant are different in the long- and short-distance parts. 
However, this effect, related to logarithms of n, can be treated consistently only in 
the context of a full renormalization group treatment. This has been done in the 
present work (see [0, |[), but not in [|I5], |T6[. As a consequence there is no analogue 
of /ifac in the present approach, while the role of /ihard in HHj is taken by the 
starting scale for the renormalization group evolution. As mentioned earlier, this 
dependence is negligible in our representation of the moments. 

(d) Hoang and this work expand the bound-state (5-functions for reasons discussed 
earlier; MY keep them unexpanded. This increases the theoretical moments signif- 
icantly at NNLO. Not expanding the bound state (5-functions would increase the 
value quoted in (|^) and (||) by almost 100 MeV. 



'As the analysis in [g^ supersedes [|l4|, we do not discuss ||lj| in detail. 
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(e) MY and Hoang use a 2-loop formula to obtain nib from thefr intermediate mass. As 
explained above, a NNLO analysis of the sum rule determines the PS (or related) 
masses with mf,af accuracy and to fully exploit this accuracy, the 4-loop relation 
between the PS and MS mass should be used. We estimated that the 3- and 4-loop 
terms decrease m?, by 45 MeV and this additional shift has been incorporated in 
(^. Employing a less accurate relation entails a corresponding loss in parametric 
accuracy of mj,, although this is a numerically small effect, if our estimate is correct. 

As in this work, MY obtain their result from an analysis of single moments (check- 
ing consistency between a set of moments), although larger moments n = 14-18 are 
used, which could be considered problematic. They use the so-called kinetic mass as 
intermediate mass definition. The analogue of 6m2{fif) in (^), needed for a NNLO sum 
rule analysis, is not yet known in this scheme; MY estimate it in the large-/3o limit, an 
additional assumption we had to use only for the 4-loop term when relating (|) to @, 
but not to extract the PS mass in the first place. MY vary the renormalization scale 
from 2 GeV to 4.5 GeV, while we would argue that, for the high moments used in [ p!5[] , 
the scale should be chosen lower. If we repeat our analysis with the same assumptions 
as those of MY, we reproduce their error estimate for the kinetic mass, which is smaller 
than the more conservative procedure that leads to (^. 

Hoang uses an analysis and error estimate that is different from the single-moment 
analysis performed by MY and in this work. Simplifying somewhat, Hoang fits the quark 
mass from the linear combination O.I2M4— 0.56M6+0.76M8— O.SlMio of moments, where 
the coefficients are determined by the covariance matrix of the experimental input data 
such as the measured leptonic widths of the six T resonances. This linear combination 
(which entails a cancellation of one part in 4000) turns out to be very insensitive to the 
renormalization scale n, yet retaining a large sensitivity to nib- Hoang then scans the 
theoretical parameter space and finds an error of only ±30 MeV for the so-called IS-mass, 
compared to the error of (^. We have repeated our analysis for this linear combination 
and obtain rrib ps(2 GeV) = 4.58 ± 0.02 GeV in this way, where the quoted error is due 
to variation of the renormalization scale only. The result is consistent with (^), but the 
error is much smaller. The central value is 50 MeV higher than the value reported in 



Sect. 6 of [|16|. Such differences can be explained by different implementations of the 
NNLO result, as discussed above. 

Several circumstances make us suspicious that the theoretical error is underestimated 
by Hoang's procedure. For example, if we increase the error of the leptonic width of the 
T(2S) by a factor of 10, or if we increase the error on the measured mass of the T(1S) 
to 5 MeV, which is still small compared to the expected theoretical error, the procedure 
chooses a linear combination that exhibits less stability in /i, and has no or two solutions 
for TTib ps(2 GeV) for some ranges of /i, even though our experimental knowledge of the 
width or mass should have no bearing on the theoretical error estimate. This remark 
may not be considered as a serious objection, because we could abandon the way the 
linear combination is chosen in |T6| and optimize it deliberately. However, even for the 



original linear combination, there is a second solution mb,ps(2GeV) = 4.86 ± 0.04 GeV 
in addition to mbPs(2GeV) = 4.58 ± 0.02 GeV, because the linear combination is no 
longer a monotonic function of rUb. The criterium of renormalization scale stability does 
not exclude obtaining solutions that differ by more than the error estimated from the 
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/i-dependence. The problem is compounded by the observation that the stabihty under 
variations of the renormahzation scale, and hence the small error obtained by Hoang, 
crucially depends on the assumption that the four moments are combined at the same 
value of the renormahzation scale fi. This is a serious assumption, in particular as the 
natural scale of the moments is 2mb/y/n. If we combine the moments at their natural 
rather than at equal scales, the stability is lost. To be fair, we should mention that 
Hoang's analysis is more involved than analysing a single linear combination, although 
the covariance matrix is such that it does indeed give most weight to a single one. 
Nevertheless, we think that the simplified discussion above emphasizes the problem with 
estimating a theoretical error in the way done in []16| . 

The final results by MY {nib = 4.2 ± 0.1 GeV) and by Hoang (mf, = 4.20 ± 0.06 GeV) 
agree with (P) within the quoted errors. However, if the 45 MeV shift were applied to 
those results, there is actually a discrepancy of about 100 MeV in the central value. 
This could be a consequence of the different representations of the moments as dis- 
cussed above. However, we find it difficult to reconcile a Wib significantly smaller than 



4.25 GeV with the analysis of Mx{is) [see ([T^) ], unless there is indeed a large positive 
non-perturbative contribution to Mx{is). 

Summary. We determined the bottom quark mass in the MS scheme and the potential 



subtraction (PS) scheme [10| at next-to-next-to-leading order from sum rules for the bb 
cross section and the mass of the T(1S) state. The results are in excellent agreement with 
each other as summarized by (|), (|^), (|l5) and (^2[). There is no systematic procedure to 
combine the two results. On the one hand, the two determinations are not independent, 
because some theoretical input is common to both. On the other hand, the dominant 
source of theoretical error is different. We therefore combine the two determinations to 
yield the PS mass 

mbps(2 GeV) = 4.59 ± 0.08 GeV (13) 

and the MS mass (at the scale of the MS mass) 

rribimb) = 4.25 ± 0.08 GeV. (14) 

The calculations that go into these results imply partial resummations of the QCD 
perturbative expansion to all orders. An important point is that a NNLO resummation 
allows us to determine the quark masses with a parametric accuracy of order rriba'^, 
i.e. the residual error scales formally as Smb/rub ~ af. In the case of the MS mass 
this requires that one controls the four-loop relation to the PS mass. We estimated the 
4-loop term, which is not yet known exactly and found that it should be very small. 

Unfortunately, the sum rule analysis yields a much less precise determination of the 
bottom quark mass than what might have been expected with NNLO accuracy. We 
identify as the reason for this the bad behaviour of the perturbative expansion for the 
leptonic width of the T resonances. The same is true for the T(1S) mass in the on-shell 
scheme. However, in this case it is understood that the large coefficients are unphysical 
and can be removed by a suitable mass subtraction procedure. If a similar mechanism 
underlied the expansion for the leptonic width, the error on the bottom quark mass could 
be reduced. In the absence of any understanding of this point, we have adopted a more 
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conservative error estimate than in previous works |]13[ |14[ |15], |16|, mainly because of a 
more generous variation of the renormahzation scale. Eq. ( |I^ ) is in agreement within 
errors, but larger than the quark mass values quoted there, but is about 200 MeV smaller 
than the one in ^^. We argued that the result of |T^ should be corrected for the large 



3-loop term in the relation between the pole mass and the MS mass. 

Eq. (0) is also in good agreement with rnb{rnb) = 4.26 ± 0.07 GeV found in [ p^ . 



This work uses the B meson mass, a lattice calculation of the (properly defined) binding 
energy of the B meson in the unquenched, two-flavour approximation to heavy quark 
effective theory, and a two-loop perturbative matching to the MS scheme. To our knowl- 
edge, this is the only other NNLO determination of the MS mass besides the sum rule 
calculations mentioned above (which, in fact, are N^LO as far as nib is concerned). How- 
ever, because of the heavy quark limit, there are l/rrib corrections, which remain to 
be estimated. Finally, the result is also in agreement with earlier, parametrically less 



accurate determinations, as for example in |29 
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